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I. INTRODUCTION 

The investigation of long-lived heavy quarkonium Be 
composed of quarks with different flavors can produce a 
significant progress in the study of heavy quark dynam- 
ics, since the variation of bound state conditions for the 
heavy quarks in various systems such as the heavy-light 
hadrons or doubly heavy mesons and baryons"'^ provides 
us with the different conditions in both the binding of 
quarks by the strong interactions and the electroweak 
decays. In addition to the rich fields of study such as the 
spectroscopy, production mechanism and lifetime there 
is a possibility to get the model-independent information 
on the CP-violating parameters in the heavy quark sec- 
tor^ 0. 

The first experimental observation of the Be meson by 
the CDF collaboration 4] confirmed the theoretical pre- 
dictions on its mass, production rate and lifetime So, 
we could expect, that an essential increase of statistics 
in the nearest future will provide us with a new battle 
field in the study of long-lived doubly heavy hadrons at 
Tevatron and LHC 0. 

The special role of heavy quarks in QCD is caused by 
the small ratio of two physical scales: the heavy quark 
mass toq is much greater than the energetic scale of the 
quark confinement Aqcq. This fact opens a way to de- 
velop two powerful tools in the description of strong inter- 
actions with the heavy quarks. The first tool is the per- 
turbative calculations of Wilson coefficients determined 
by the hard corrections with the use of renormalization 
group improvements. The second instrument is the Op- 
erator Product Expansion (OPE) related with the small 
virtuality of heavy quark in the bound state, which re- 
veals itself in many faces such as the general expansion of 
operators in inverse powers of heavy quark mass and the 
QCD sum rules. A specific form of OPE is the application 



for the heavy quark lagrangian itself, which results in the 
effective theories of heavy quarks. The effective theory is 
constructed under the choice of its leading term appropri- 
ate for the system under study. So, the kinetic energy can 
be neglected in the heavy-light hadrons to the leading or- 
der. The corresponding effective theory is called HQET 
12 . In the doubly heavy mesons the kinetic energy is of 
the same order as the potential one, while the velocity of 
the heavy quark motion is small, and we deal with the 
nonrelativistic QCD 9] and its developments by taking 
into account the static energy in the potential NRQCD 
in a general form (pNRQCD) or under the correla- 
tion of quarkonium size and the time for the formation 
of bound system with the improved scheme of expansion 
in the heavy quark velocity (vNRQCD) [l^- The per- 
turbative QCD remains actual in the effective theories, 
since it is necessary for the calculation of Wilson coeffi- 
cients determined by hard corrections. More definitely, 
the perturbative calculations determine both the match- 
ing of Wilson coefficients in the effective theory with the 
full QCD and the anomalous dimensions resulting in the 
evolution of the coefficients with the variation of normal- 
ization point. In this respect we have to mention that the 
pNRQCD results on the static potential and the mass- 
dependent terms in the heavy quark-antiquark energy 
were confirmed by vNRQCD after appropriate limits and 
some calculational corrections in vNRQCD. In addition, 
the pNRQCD is a powerful tool in the studies of both 
the spectroscopy and the heavy quarkonium decays [lll| . 

In contrast to the Wilson coefficients, the hadronic ma- 
trix elements of operators composed by the effective fields 
of nonrelativistic heavy quarks cannot be evaluated in the 
perturbative manner. In this paper we use the potential 
model for such the calculations in the case of leptonic 
constant for the mixed-flavor heavy quarkonium Be- The 
matching procedure is performed with the two-loop ac- 
curacy available to the moment. 
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^ See the review of physics with the baryons containing two heavy 
quarks in 

^ The extraction of angle 7 in the unitary triangle derived from 
the charged current mixing in the heavy quark sector can be 
obtained from the decays of doubly heavy baryons, too, in the 
same manner la. 



II. LEPTONIC CONSTANT OF Be IN THE 
POTENTIAL APPROACH 

In the NRQCD approximation for the heavy quarks, 
the calculation of leptonic constant for the heavy quarko- 
nium with the two-loop accuracy requires the matching 
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of NRQCD currents with the currents in full QCD, 



NRQCD 



where we have introduced the following notations: Qi,2 
are the relativistic quark fields, x s-^d are the nonrel- 
ativistic spinors of anti-quark and quark, v is the four- 
velocity of heavy quarkonium, so that 



(1) 



where the scale /ihard gives the normalization point for 
the matching of NRQCD with full QCD, while /ifact de- 
notes the normalization point for the calculations in the 
perturbation theory of NRQCD. 

For the pseudoscalar heavy quarkonium composed of 
heavy quarks with the different flavors, the Wilson coef- 
ficient K, is calculated with the two-loop accuracy 
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as''(A'hard) 
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A^(/^hard; /^fact) 



and ci,2 are explicitly given in Refs. J3| and respec- 
tively. The anomalous dimension of factor /C(^hard; Mfact) 
in NRQCD is defined by 



(iln/C(^hard;M) fcC^iP') 
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whereas the two-loop calculations'^ give 
7[i] = 0, 

7[2] 
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where r denotes the ratio of heavy quark masses. The 
initial condition for the evolution of factor /C (/ihard; Mfact) 
is given by the matching of NRQCD current with full 

QCD at /Xfact = A'hard- 

The leptonic constant is defined in the following way: 



{Q\jd^^\QQ)=V,fQQMQQ. 



(6) 



In full QCD the axial vector current of quarks has 
zero anomalous dimension, while in NRQCD the current 
jjjNRQCD j^g^g ^Yie nonzero anomalous dimension, so that 
in accordance with C])-©, we find 



'QQ 



(7) 



We use ordinary notations for the invariants of SU(Nc) repre- 

N^ — l 1 
sentations; Cp = -fj^ — , Ca = -'Vc, Tp = ^, nj is a number of 

"active" light quark flavors. 



where, in terms of nonrelativistic quarks, the leptonic 
constant for the heavy quarkonium is given by the well- 
known relation with the wave function at the origin 



NRQCD 

QQ ^ 
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and the value of wave function in the leading order is 
determined by the solution of Schrodinger equation with 
the static potential, so that we isolate the scale depen- 
dence of NRQCD current in the factor -4(/i), while the 
leptonic constant Iq^'^^ is evaluated at a fixed normal- 
ization point fi — fio, which will be attributed below. It 
is evident that 

^QQ = /Oo'^^-^^^f-*) • ^(Mhard; Mfact), (9) 

and the anomalous dimension of .4(/ifact) should compen- 
sate the anomalous dimension of factor IC (/ihard ; /ifact ) , so 
that in two loops we have got 



d\r^A{^J,) 
dln/i 



-7[2] 



47r 



(10) 



The physical meaning of A{p) is clearly determined by 
the relations of lO and l^: this factor gives the nor- 
malization of matrix element for the current of nonrel- 
ativistic quarks expressed in terms of wave function for 
the two-particle quark state (in the leading order of in- 
verse heavy quark mass in NRQCD). Certainly, in this 
approach the current of nonrelativistic quarks is factor- 
ized from the quark-gluon sea, which is a necessary at- 
tribute of hadronic state, so that, in general, this phys- 
ical state can be only approximately represented as the 
two-quark bound state. In the consideration of leptonic 
constants in the framework of NRQCD, this approxima- 
tion requires to introduce the normalization factor A{fJ.) 
depending on the scale. 

The renormalization group equation of (|10|l is simply 
integrated out, so that 



A{^l) = Ai^lo) 



/3o + /3i 



An 



f3o + Pi 



'(a^o) 



(11) 



where /3o = ^CA-jTpnf, and /3i = ^-Cl- AC pTpn f - 
^CaTpUj. a constant of integration could be defined 
so that at a scale fj,Q we would get -4(/io) — 1. Thus, in 
the framework of NRQCD we have got the parametric 
dependence of leptonic constant estimates on the scale 
fiQ, which has the following simple interpretation: the 
normalization of matrix element for the current of non- 
relativistic quarks at fiQ is completely given by the wave 
function of two-quark bound state. At other /z ^ /^o 
we have to introduce the factor A{fJ,) 7^ 1, so that the 
approximation of hadronic state by the two-quark wave- 
function becomes inexact. 
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Following the method described in 0,^^, we estimate 
the wave function of be quarkonium in the nonrelativistic 
model with the static potential given by 0|, so that 

V^\^{0)l^\ 1.26 GeV^/2. (12) 

In a wide region of distances between the quarks p, the 
static potential of free of infrared singularity can be 
matched with the two-loop perturbative potential V^crt 
calculated in pJt QM ■ ^^^^ '^^y difference between 
the potentials 

SV{Hsolt) = Vkko(p) - Vport(/3; ^J■solt) (13) 

depends on the normalization point fieoft in the pertur- 
bative QCD of static sources. It can be numerically ap- 
proximated by the expression 

^^(A^soft) ^ 0.99 
IGeV Ai,oft/(lGeV)-0.5' 

reflecting the infrared uncertainties in the perturbative 
calculations. 
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FIG. 1: The leptonic constant of ground pseudoscalar state 
in the system of heavy quarkonium be is presented versus the 
soft scale of normalization. The shaded region restricted by 
curves corresponds to the change of hard scale from /Xhard = 3 
GeV (the dashed curve) to Hhard ~ 2 GeV (the solid curve) 
with the initial condition for the evolution of normalization 
factor ^(pifact) posed in the form of ^(1.2 GeV) = 1 and 
^(1. GeV) = 1, respectively, in the matrix element of current 
given in the nonrelativistic representation. The horizontal 
band is the limits expected from the QCD sum rules [20| and 
scaling relations for the leptonic constants of heavy quarkonia 
PH . In the cross-point, the leptonic constant of Be weakly de- 
pends on the parameters given by the hard scale of matching 
as well as the scale of the initial normalization. 

The masses of heavy quarks used in the potential 
model of (15j 

= 1.468 GeV, = 4.873 GeV, (14) 



are consistent with the values of running masses as was 
shown in [T5l |. The cancellation of renormalon in the 
sum of the perturbative static potential and pole masses 
of heavy quarks takes place, so that'' 

TOpolG(Msoft) + ^ SV{flsoft)- 

In what follows we put the normalization point of per- 
turbative potential /isoft = Mfact- 

The result of calculation for the leptonic constant of 
Be in the potential approach is shown in Fig. We 
chose the values of /^hard in a range, so that the stable 
value of leptonic constant would be posed at the scale 
A'fact < ?7ic°'°: which is the condition of consistency for 
the NRQCD approach. Then one can observe the po- 
sition of /ifact ~ 1.5 — 1.6 GeV, where the estimate of 
/b^ is independent of the choice of /ihard and /io, and the 
maximal values at the different choices of /ihard are also 
close to each other. 

/flc, GeV 
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FIG. 2: The leptonic constant of ground pseudoscalar state 
in the system of heavy quarkonium be is presented versus the 
hard scale of normalization at /ifact = 1-5 GeV. 

We have found that in order to reach the minimal de- 
pendence of leptonic constant on the hard scale, the nor- 
malization point fiQ should correlate with /ihard at fixed 
/ifact- For illustration we show the result at /ifact = 1-5 
GeV and 

/io = 1.6 GeV - (51^ (/ihard) 

in Fig. m 

The full dependence on the scales is presented in Fig. 
13 We can see that there is a saddle stable point for the 



* The accuracy of ex pression below, i.e. a possible additive shift, 
is discussed in Il9l . 
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FIG. 3: The leptonic constant of ground pseudoscalar state 
in the system of heavy quarkonium be is presented versus the 
soft and hard scales, /ifact and ^hard- The initial condition 
for the evolution of normalization factor ^(/Jfact) in the ma- 
trix element of current in the nonrelativistic representation is 
given by ^(1.6 GeV - 5V(Mhard)) = 1. 

leptonic constant depending on two scales: the normal- 
ization point of perturbative calculation for the Wilson 
coefficient in NRQCD and the matching point of NRQCD 
with the full QCD. This stability implies that in the vicin- 
ity of saddle point the contribution of higher corrections 
in as is not significant, while at other values of scales 
these corrections should be taken into account. 
The final result of two-loop calculations is 

/b^ = 395 ± 15 MeV. (15) 



The result on /b^ is in agreement with the scalin g re - 
lation derived from the quasi- local QCD sum rules |2ll| . 
which use the regularity in the heavy quarkonium mass 
spectra, i.e. the fact that the splitting between the 
quarkonium levels after the averaging over the spins of 
heavy quarks weakly depends on the quark flavors. So, 
the scaling law for the S-wave quarkonia has the form 

/« / mi + ma y ^ c 

where n is the radial quantum number, mi^2 sue the 
masses of heavy quarks composing the quarkonium, /ii2 
is the reduced mass of quarks, and c is a dimensional con- 
stant independent of both the quark flavors and the level 
number n. The value of c is determined by the splitting 
between the 25 and IS levels, or the average kinetic en- 
ergy of heavy quarks, which is independent of the quark 
flavors and n with the accuracy accepted. The accuracy 
depends on the heavy quark masses, and it is discussed 
in |2ll | in detail. The parameter c can be extracted from 
the known leptonic constants of ip and T, so that the 
scaling relation gives 

/s. « 400 MeV 

for the vector state. The difference between the lep- 
tonic constants for the pseudoscalar and vector IS'-states 
is caused by the spin-dependent corrections, which are 
small. Numerically, we get — /bJ/Jb* < 3%, hence, 
the estimates obtained from the potential model and the 
scaling relation is in a good agreement with each other. 



It should be compared with the estimate of potential 
model itself without the matching 

which indicates the magnitude of the correction about 
20%. Furthermore, the calculations in the same potential 
model with the one-loop matching gave 

where the uncertainty is signiflcantly greater than in the 
two-loop procedure, since at the one-loop level we have 
no stable point in the scale dependence of the result. 
Therefore, in contrast to the discussion given in we 
see that the correction is not crucially large, but it is 
under control in the system of Be- The reason for such 
the claim on the reliability of result is caused by two cir- 
cumstances. First, the one-loop anomalous dimension of 
NRQCD current is equal to zero. Therefore, we start the 
summation of large logs in the framework of renormaliza- 
tion group (RG) with the expressions in (|10|l and (|ll|l . 
Second, after such the summation of large logs the three- 
loop corrections could be considered as small beyond the 
leading RG logs. 



III. CONCLUSION 

In this paper we have calculated the leptonic constant 
for the ground state of be system in the framework of the 
QCD-motivated potential model used for the estimate of 
nonperturbative hadronic matrix element for the effective 
fields of nonrelativistic quarks. The Wilson coefficient 
calculated with the two-loop accuracy ,,14] has been im- 
plemented to relate the currents of heavy quarks in the 
full QCD and NRQCD. It depends on two scales: the 
matching point and the normalization scale in NRQCD. 
The stable point in the scale parameter space has been 
observed, which makes the estimate to be reliable. We 
have found that the two-loop corrections are under con- 
trol in the Be meson. 

The numerical value of /b^ is in agreement with the 
estimates obtained in the framework of QCD sum rules 
[23, potential models ^ and from the scaling relation 
|2ll |. The expected systematic error coming from the us- 
age of potential model is reduced due to the QCD motiva- 
tion consistent with the measurements of QCD coupling 
constant at large virtualities, the asymptotic freedom up 
to three loops and the fitting the mass spectra of char- 
monia and bottominia, so that the systematic accuracy 
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is of the same order as the uncertainty related with the 
variation of parametric scales shown in l|15|l . 
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